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THE STRUCTURAL EVOLUTION OF MILKY WAY-LIKE STAR FORMING GALAXIES SINCE 
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ABSTRACT 

We follow the struetural evolution of star forming galaxies (SFGs) like the Milky Way by selecting 
progenitors to z ~ 1.3 based on the stellar mass growth inferred from the evolution of the star 
forming sequence. We select our sample from the 3D-HST survey, which utilizes spectroscopy from 
the HST WFC3 G141 near-IR grism and enables precise redshift measurements for our sample of 
SFGs. Structural properties are obtained from Sersic profile fits to CANDELS WFC3 imaging. The 
progenitors of z = SFGs with stellar mass M = 10^"'^ Mq are ^ 2 times less massive at z ~ 1. This 
late-time stellar mass assembly is consistent with recent studies that employ abundance matching 
techniques. The descendant SFGs at z ~ have grown in half-light radius by a factor of ~ 1.4 
since z ^ 1. The half-light radius grows with stellar mass as re oc M"-^^. While most of the stellar 
mass is clearly assembling at large radii, the mass surface density profiles reveal ongoing mass growth 
also in the central regions where bulges and pseudobulges are common features in present day late- 
type galaxies. Some portion of this growth in the central regions is due to star formation as recent 
observations of Ha maps for SFGs at z ~ 1 are found to be extended but centrally peaked. Connecting 
our lookback study with galactic archeology, we find the stellar mass surface density at i? = 8 kpc to 
have increased by a factor of ~ 2 since z ~ 1, in good agreement with measurements derived for the 
solar neighborhood of the Milky Way. 
Subject headings: galaxies: structure — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

The formation of late-type star forming galaxies 
(SFGs) like the Milky Way has been studied through sev- 
eral complementary approaches in order to understand 
the assembly of the disk and bulge. On the observa- 
tional front, detailed studies of the ages, metallicities, 
and kinematics of stellar populations within the Milky 
Way itself afford a unique vantage point for viewing the 
assembly of such late- type galaxies. This extremely infor- 
mative approach is limited in completeness and sample 
size and must contend with a reconstruction of the global 
star formation and struc tural assembly history, which is a 
comp lex task (see, e.g., iRoskar et al.ll2008t iRix fc Bovvl 
120131 and references therein). On the theoretical front, 
the formation of Milky Way-like SFGs has posed a chal- 
lenge for the numerous simulations that aim to produce 
a structurally realistic analog at z ~ with a disk and 
an em bedded bulge or pseudobulge (jScannapieco et ahl 
I2012D . In particular, this latter work shows how many 
of the simulations form massive bulges at early times, 
lack prominent disks at late times, and produce widely 
varying physical properties for the descendant galaxy 



Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science In- 
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26555. 



at z = 0. This largely refiects the challenge of cor- 
rectly modeling the "sub-grid" physical processes that 
impact the baryons. Recently, however, there has been 
some success among a few simulations in producing late- 
type SFGs like the Mil k Way within sp ecific halos (e.g., 
Guedes et aljl'2011; S tinson et al.ll2013l ). 

In this Letter, we employ a lookback study that com- 
plements the galactic archeology and theoretical methods 
above in order to understand how late- type SFGs like the 
Milky Way were assembled. We select star forming pro- 
genitors using the stellar mass growth inferred from the 
evolution of the star forming sequence and analyze their 
structural evolution using HST imaging. This method 
for computing th e mass growth of SFGs has been dis- 
cussed in detail bv lLeitne"ir(|2012| ) and complements other 
studies which conn ect progenitors and descendants using 
numb er densities (jvan Dokkum et ahl 120101 iPatel et al.l 
l2013f l . Given the distinct formation history of SFGs from 
quiescent galaxies (QGs), we use the former method here 
as it directly traces the sta r forming progenito r s of ga lax- 
ies like the Milky Way. In Ivan Dokkum et ahl (|2013f ) , we 
use the number density approach to study progenitors of 
galaxies of all types with the mass of the Milky Way. 

We assume a cosmology with Hq = 70 km s~^ Mpc^^, 
fiM = 0.30, and f^A = O-'^'Q- Ste llar masses are based on 
a Chabrier IMF (jChabrieill2003f ). Magnitudes are in the 
AB system. 
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Figure 1. (a) Evolution of SFGs in tiie SSFR-mass plane witli final mass at 2 ~ of _A/ = 10^"'^ Mq as computed by ILeitnerl II2012I ) 
using tile SSFFl-mass relations at different redshifts (dashed lines) from Karim et al. (2011). At a given redshift, new stellar mass is added 
according to tiie position of the SFG on the star forming sequence and mas s loss from stellar evolution is also accounted for. The slight 
offset between the extrapolated 2 = relation from that of lSalim et al.l l|2007l ^ does not impact our results (see text). The Milky Way (gray 
circle) lies within the observational scatter of the 2 = relation, (b) SFH of the SFGs tracked in panel (a). The data points with error 
bars represent IR+ UV SFRs for a subset of our sample with deep MIPS: there is good general agreement with the radio based SFRs from 
IKarim et al.l 120111 ). (c) Mass growth history of the SFGs tracked in panel (a). We select galaxies at the indicated masses and redshifts in 
this Letter. Since z = 1, the SFGs grow in ste llar mass by a factor of ~ 2. For comparison, the SFH and mass growth for galaxies with 
the same final mass from lBehroozi et al.l II2012I) is shown by the dashed red li ne. There is excellent agreement between the observationally 
motivated method used here and the more theoretically motivated method of lBehroozi et aLl 120121 ). 



2. DATA AND ANALYSIS 

We employ data from the 3D-HST survey (w2.1) to 
carry out our analysis. The observations and data reduc- 
tion p rocedures are explained in detail in lBrammer et aLl 
(|20T1 and Skelton et al. (in prep.). The HST WFC3 
G141 near-IR grism observations are the centerpiece of 
3D-HST and cover the CANDELS fields. In this work, we 
use the three fields with currently available WFC3 based 
structural parameters from the literature: COSMOS, 
GOODS-S, and UDS. Grism reds hifts were measured us - 
ing a modified version of EAZY (jBrammer et al.ll2008[ ). 
The procedure fits template SEDs to both the photom- 
etry and grism spectroscopy, enabling precise redshift 
measurements for our sample of SFGs at hi gh redshift. 

Ste llar masses were measured with FAST (jKriek et al.l 
l2009f l . We determine the stellar mass limit at our highest 
redshifts {z = 1.4) to be M ~ 10^® M^;, (see Figure [He)) 
using a similar technique to that of iMarchesini et al.l 
(|2009f ). This conservative limit accounts for galaxies with 
high M/L such as reddened SFGs. Rest-frame U — V and 
V — J colors were measured with EAZY in order to dis- 
tinguish SFGs from QGs (Section [Sj). 

Structu ral properti es were obtai ned from 
van d or Wei et all (poTl who used GALFIT (|Peng et al.l 
200l to fit Sersic profiles to t he CANDELS ffgT 
WFC3 Ji25 and JJia o imaging ()Grogin et al.l 120111: 
iKoekemoer et al.1 120111) . The effective radii reported 
here are circularized. The J125 imaging is used at 
0.25 < z < 1 while the iJieo imaging is used at 
1 < z < 1.4, with our median galaxy sampling rest- 
frame - 7100 A. The PSF FWHM/2 of the iJieo 
imaging is ^ 0'.'09, corresponding to a physical radius at 
our highest redshifts of ~ 0.8 kpc. 

We employ an SDSS reference sample at z = 0.05 using 
stellar masses and SFRs from the MPA/JHU catalogs 
(jBrinchmann et al.li2004') and z-band Se rsic fits from the 
NYU-VAGC (IBlanton et aLll2003Ll2005l) . 



3. SELECTION 

3.1. Stellar Mass Growth History from an Evolving Star 
Forming Sequence 

In this work, we trace the formation of SFGs with a fi- 
nal mass slightly below that of the Milky Way (~ 0.2 dex) 
so that we can compare various properties to galaxies 
of the same final mass from other works that employ 
different progenitor -descendant linking methods (e.g., 
iBehroozi et aL|[2012l ) . We note however that we arrive at 
similar qualitative conclusions when tracing progenitors 
of Milky Way-mass SFGs. 

Selecting progenitors of z ~ SFGs requires knowl- 
edge of their mass growth history so that one can se- 
lect galaxies of the proper progenitor mass at a given 
redshift. For galaxies that assemble most of their stars 
from in-situ star formation, one can infer the mass 
growth from the evolution in the SFR-mass relation (or 
SSFR-mass). The Milky Way, with a stellar mass of 
Af ~ 5 X 10^° Mq CHammcr et al."2007^ and a SFR of 
1.9±0.4 Mq yr-i (Chomiuk & Povich 2011), faUs within 
the observa t ional scatter of the SSFR-mass relation of 
iSalim et"aLl (j2007l their Eq. 11, corrected for evolution 
to z = assuming Equation [1] here; see Figure jlfa)). 
This star-forming sequence has been extens ively studied 
to high redshifts (e.g., [N ocskc ct al. 20 071: IKariin et all 
120111: IWhitaker et al.|[2Q'l2l: IFumagalh et aLll2012D . The 
method for determining the mass growth from this re- 
lation is simple: for a given redshift interval, new stel- 
lar mass is added to the existing mass based on the lo- 
cation of a galaxy in the SSFR-mass plane and mass 
loss from stellar evolution is accounted for from simple 
stellar population models. The assumption that most 
nearby SFGs were star forming at high redshift is sup- 
ported by the small scatter about the SSFR-mass rela- 
tion {a ^ .17 dex after accounti ng for observational un- 
certainties; |Whitater]eOT][20T^ as well as by the vastly 
differ ent structural pro perties between SFGs and QGs 
(e.g., iFranxet al.ll2008f ). 
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Figure 2. (a)-(e) Rest-frame U — V versu s V — J at different redsfiifts, with galaxies selected to have the progenitor mass (±0.1 dex) 
from Figure[llc). The Williams et al.l 120091) boundary separating QGs from SFGs (bottom right) is used to select the latter. These SFGs 
are color-coded according to their Sersic indices, (f) WFC3 G141 median stacked spectra of SFGs for redshift bins with observable Ha. 

In th i s wor k, we use the mass growth computed by 
iLeitneii (I2012D. who derived it from the SSFR-mass re- 
lations of iKarim et"al] ()20lH) . Figure [Ha) shows SSFR- 
mass rel ations at different re dshifts (dashed lines) mea- 
sured by IKarim et al.l (|2011) from 1.4 GHz radio stack- 
ing. They employ the parameterization: 



SSFR ozM'^{l + zy 



(1) 



where (3 — —0.35 and n = 3.45. Though the z = 
relation based on E quation [T] is an extrapolation of the 
IKarim et al.l ()2011[ ) data, it s slight difference from lower 
redshift observations (e.g.. ISalim et al]|2007D does not 
significantly impact the mass growth at these late times. 
The track shown in Figure [Ija) (black curve) indicates 
the trajectory in the SSFR-mass plane for the SFGs se- 
lected for study in this work with final mass at z = 
of M = 10^°-^ Mq. The corresponding SFH is shown 
in panel (b) and the mass growth with redshift in panel 
(c). The points with error bars (including a 30% sys- 
tematic uncertainty) in panel (b) indicate the median 
SFR (IR+UV) of our sample in the GOODS fields (i.e., 
deep MIPS) for our C/FJ-selected (see below) SFGs at 
the corresponding redshifts and masses from panel (c). 
While SFR estimates can v ary due to s ystematics be- 
tween different tracers (e.g.. iWuvts et al..,2011; iLeitner. 



120121 ) there is good general a greement between o ur SFR 
measurements and those of IKarim et al.l ()201lD . The 
dashed red curves in panels (b) and (c) show the cor- 
responding result for the same final mass galaxy from 
iBehroozi et al.l (2012) whose results were derived from a 
more theoretical approach that combined merger trees 
from dark matter simulations with observational con- 
straints. The latter is consistent with our observationally 
motivated method. 

The results above suggest significant stellar mass as- 
sembly in our SFGs below z < 1 (factor of ~ 2.2). 
Late assembly at z < 1 is also found by othe r abun- 
dance matching works (e.g.. iMoster et al1l2013f ). SFHs 
of the Milky Way disk derived from stellar proper- 
ties also suggest signifi cant late time ass embly (e.g., 
iRocha-Pinto et al.l I2OOOI : lAumer fc Bi"^!^ [2001 . Fi- 
nally, hydrodynamical simulations also appear to favor 
late time assembly for Milky Way-like galaxies given that 
earlier formation fails to p roduce a prominent stellar disk 
([Scannapieco et aL[|2012t ). 

3.2. Selection of Star Forming Progenitors to z '^ 1.3 

Given the observational scatter about the SSFR-Mass 
relation, in practice we select all SFGs in a given red- 
shift and mass bin. Figure [2] shows rest-frame U — V 
vs. V — J diagrams for galaxies in narrow mass bins 
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Figure 3. Example SDSS i-band and HST WFC3 postage stamps for progenitors of SFGs with a final mass at 2 = of Af = 10^*^'^ Mq. 
Each stamp is ~ 30 kpc on a side. The median stellar mass decreases to high redshift according to Figure [TJc) where the sizes appear 
smaller. 



(±0.1 dex) centered on the progenitor mass at a given 
redshift. We use this UVJ selection to identify SFGs, 
which occupy the bottom right portion o f these diagrams 
(iWihiams et aLl 120091: Fatel et all |20T2| ). The descen- 
dants, which are more massive, become redder toward 
low redshift likely due to their aging stellar populations. 
In the SDSS sample, galaxies with SSFR > 10"" yr^^ 
were selected as SFGs. Figure [3] shows example SDSS 
i-band and HST WFC3 postage stamps for progenitor 
SFGs at different redshifts. 

4. RESULTS 

4.1. Structural Evolution 

We examine the evolution in various structural param- 
eters for our sample of progenitor SFGs in Figure |4l 
Panel (a) shows that most of the progenitors at high 



redshift have close to exponential profiles while those at 
lower redshift have slightly higher Sersic indices. The me- 
dian axis ratio in panel (b) remains relatively constant 
at b/a w 0.60, a low value which for a population of ran- 
domly inclined disks would imply an intrinsic axis ratio 
of ~ 0.2. The low Sersic indices at higher redshifts and 
low axis ratios therefore point to the progenitor SFGs as 
being disks. 

Half-light radii provide a first order view into the dis- 
tribution of stellar mass for galaxies in our sample. Fig- 
ure \M,c) shows the evolution of the median half-light ra- 
dius with redshift. Though we follow the median of a 
given property, it is important to note that at a given 
redshift and stellar mass, SFGs display a diversity of 
property values. The black dashed line represents a fit 
to the data of the form: 
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Figure 4. (a) Median Sersic index versus redshift for progenitor SFGs. (b) Median axis ratio versus redshift. (c) Median half-light radius 
versus redshift. The error bars represent the bootstrapped uncertainty on the median. Fits of the form (1 + z)" are indicated for the Sersic 
indices and half-light radii, both of which increase toward low redshift. (d) Evolut i on of re in the size-mass plane. The dotted blue line is 
the z ~ SDSS size- mass relation (z-band) for late- typo galaxies from lShen et al] II2003I) . The dashed line is a fit to the data of the form 
re tx M", where a = 0.29 ± 0.08. Since z = 1, the half-light radii of the SFGs have grown by a factor of ~ 1.4. 

where a — 0.29 ± 0.08. Interestingly, the value of a mea- 
sured here for SFGs is much sm aller tha n that fo r QGs 
re=l3{l + z)°' (2) (ann » 2: l^anDokkumeFan [2010: Pat el et all [2013V 

SFGs generally appear to evolve near the local scaling 
relation below z < 1.3, consistent with slow evolution in 
the zeropoint of the size-mass relation for SFGs (van der 
Wei et al., in prep). 



where /3 = 3.5 ± 0.3 kpc and a = -0.48 ± 0.15. We note 
that all of the SDSS data points in Figure S] include a 
10% error added in quadrature to account for systematics 
between different measurement methods. Since z = 1, 
the median half-light radius for these SFGs has grown 
by a factor of ^ 1.4, indicative of some amount of inside- 
out growth. 

Figure lUJd) shows the evolution of the progenitors in 
the size-mass plane. For reference, the SDSS z ~ 



size-mass relation for late- type galaxies from iShen et al.l 
([2003) is indicated by the dotted blue line. The dashed 
line represents a fit to our sample of the form: 



re oc M" 



(3) 



4.2. Stellar Mass Growth in the Central and Outer 
Regions at z < 1.3 

Stellar mass surface density profiles provide a detailed 
look at the distribution of mass within galaxies. Fig- 
ure [5ja) shows the median combined mass surface den- 
sity profiles for our sample of SFG progenitors. These 
profiles wer e cons tructed in a similar manner to that of 
iPatel et all (|2013[ ). where the best-fit single component 
Sersic profiles were stacked to create the median light 
profile. A single component Sersic profile is generally 
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Figure 5. (a) Median combined mass surface density profiles. The dotted portion of each curve indicates where the uncertainty becomes 
larger than > 20%. The shaded gray region indicat es the maxi mum WFC3 PSF FWHM/2. The solar symbol indicates the mass surface 
density in the solar neighborhood as computed by IBovv et al.l (2012) but scaled down by 0.2 dex to account for the difference in mass 
between the Milky Way and our 2: ~ descendants, (b) The ratio of the mass profiles at higher redshifts to the SDSS mass profile. Shaded 
region indicates the uncertaint y. The factor of 2.4 ili 0.3 mass growth from z ~ 0.9 to 2; ~ at _Ro = 8 kpc agrees well with measurements 
in the solar neighborhood from lAumer fc BinnevI J 2009) . (c) Median stellar mass growth for the central (R <1 kpc, red) and outer regions 
(-R > 2 kpc, blue). The total stellar mass growth is shown in black. While most of the new mass growth since 2 ~ 1 has taken place in the 
outer regions, mass has also been added in the central regions of these SFGs where bulges and pseudobulges are common feature s. Some 
porti on of the growth in the central regions is due to star formation as Ha maps for SFGs at 2 ~ 1 are centrally peaked (Nelson et al.l 

mm. 



found to be a good representation of t he light profile 
for individual galaxies at high redshifts (jSzomoru et al.1 
120121: iMosleh et al.llMl) . Prior to stacking, each Hght 
profile was converted into a mass profile by normalizing 
the light within i? = 20 kpc to the total stellar mass 
of eac h galaxy. Whil e this method ignores MjL gra- 
dients, iSzomoru et al.l (|2013l ) find such gradients in the 
rest-frame g-band to be small for higher redshift galax- 
ies and that they do not evolve (here, we trace a redder 
rest-frame i-band). We plan to explore the full impact of 
MjL gradients in foUowup work. Figure [Sjb) shows the 
ratio of the mass profiles at higher redshifts to the SDSS 
mass profile. Below 2: < 1.3, stellar mass is continually 
built up at all radii. 

Figure[5l[c) shows the growth in the projected mass in- 
side (red) and outside (blue) of i? = 2 kpc (i.e., w Tg at 
z ^ 1.3). Clearly more mass has assembled at larger radii 
since z '^ 1.3, leading to larger re toward low redshift. 
However, the amount of stellar mass in the central re- 
gions, where bulges and pseudobulges are co mmon struc- 
tural features in nearby late-type galaxies ( W einzirl et al.l 
I2009f ). has also increased. Some portion of the growth in 
the central regions is due to star formation as Ha maps 
for S FGs at z '^ 1 are centrally peaked (jNelson et al.l 
l2013f ). This is qualitatively consistent with recent ob- 
servations of stars in the Milky Way bulge that display 
a wide range of metallicities and ages, implying an ex- 
tended formation history (Bensbv et a l. 2013). The on- 
going mass assembly i n the central regions may point to 
secular processes (e.g.. lKormendv fc Kennicutigl2004f ) as 
an important channel for bulge growth for SFGs in the 
stellar mass regime studied here. 



5. DISCUSSION 

5.1. A Comparison with the Stellar Mass Growth in the 
Solar Neighborhood 

Given the significant mass assembly found at large 
radii, we compare our results in such a region of our 
galaxy that has been well documented, the solar neigh- 
borhood. Though we caution that the Milky Way is just 
one such S FG, and may not be one that is an archetypal 
late- type (jHammer et al.|[2007f ). this analysis neverthe- 
less provides an intriguing comparison between our look- 
back study and galactic archeology. The solar symbol 
in Figure [51(a) indicates the stellar mass surface d ensity 
at th e solar radius {Rq — 8.0 kpc; [V anhoUcbcke et al.l 
|2009[ ) for the Milky Way from iBovv et al.. (,201^ but 
scaled down by 0.2 dex to account for the difference in 
mass with our descendant SFGs at z '^ 0. The sys- 
tematic uncertainty in this correction may account for 
the slight offset from the SDSS mass profile, though the 
rough agreement is still remarkable. lAumer fc BinnevI 
(|2009f ) estimate a mean formation time for stars in the 
solar neighborhood that corresponds to z ~ 0.9. Assum- 
ing this mean is close to the median formation redshift, 
therefore implying a factor of ^ 2 growth below z < 0.9, 
this estimate is close to the factor of ~ 2.4 ± 0.3 growth 
in our mass surface density from z ~ 0.9 (green line) to 
the SDSS mass profile at i? = 8 kpc. 

5.2. Caveats 

Two caveats to the analysis presented here warrant 
some consideration. First, the contribution of stars 
formed ex situ to the stellar mass growth is an uncer- 
tain quantity, though there is ev idence that it i s mili imal. 
Both lBehroozi et all (j2012[ ) and lMoster et all ()2013l ) find 
that for local halos of mass Mh ~ 10^^ Mq (hosting cen- 



trals of stellar mass M = lO^"'^ Af©), lit tle stellar mass 
was a ccreted. At the lowest redshifts, iBehroozi et alj 
(|2012f) find a growing contribution, but this is likely 
driven by the dominance of QGs in their halos at low red- 
shift, which primarily grow from mergers. Major mergers 
are likely rare in our sample at z < 1 as such events tend 
to destroy disks, contrary to what is observed at low red- 
shift (Figure [S]). 

Second, the QG fraction increases toward low redshift 
and as a result not all SFGs at 2: ~ 1.3 will be progenitors 
of SFGs at low redshift as some will have quenched. This 
would impact our results if there is a relation between the 
structure of high redshift SFG progenitor candidates and 
low redshift QG descendants, as the former would bias 
our measurements for the median structural property. 
Exploring these caveats is beyond the scope of this Letter 
but will be addressed in future work. 

5.3. Conclusions 

In this Letter, we have used the evolution in the SSFR- 
mass relation to determine the stellar mass growth (e.g., 
ILeitneri [20TI) for progenitors of SFGs like the Milky 
Way. We used 3D-HST redshifts and photometry to se- 
lect SFGs of the appropriate progenitor mass at different 
redshifts back to z - 1 and HST CANDELS imaging to 
follow their structural evolution. As these SFGs grew in 
stellar mass by a factor of ~ 2 since z ~ 1, most of the 
new stellar mass assembled in the outer regions, lead- 
ing to a factor of '--^ 1.4 increase in the half-light radius. 
The mass surface density profiles indicate ongoing stellar 
mass assembly also in the central regions where bulges 
and pseudobulges are common features in spirals. We 
also found good agreement between the mass growth at 
R = 8 kpc in our lookback study and that for the solar 
neighborhood of the Milky Way. 
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